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Constant-pressure molecular-dynamics simulations of phospholipid membranes in the fluid La 
phase reveal strong correlations between equilibrium fluctuations of volume and energy on the 
nanosecond time-scale. The existence of strong volume-energy correlations was previously deduced 
indirectly by Heimburg from experiments focusing on the phase transition between the La and 
the L/s phases. The correlations, which are reported here for three different membranes (DMPC, 
DMPS-Na, and DMPSH), have volume-energy correlation coefficients ranging from 0.81 to 0.89. 
The DMPC membrane was studied at two temperatures showing that the correlation coefficient 
increases as the phase transition is approached. 



Biological membranes are essential parts of living cells. 

I They not only act as passive barriers between outside and 
inside, but also play an active role in various biological 
mechanisms. The major constituent of biological mem- 
branes are phospholipids. Pure phospholipid membranes 
often serve as a models for the more complex biological 
membranes. Close to physiological temperatures mem- 
branes undergo a transition from the high-temperature 
fluid La phase (often referred to as the "biologically rele- 

] vant phase") to a low-temperature ordered gel phase Lp. 
In the melting regime response functions such as heat ca- 
pacity, volume-expansion coefflcient, and area-expansion 
coefficient increase dramatically. Also, the characteris- 

I tic time for the collective degrees of freedom increases 
and becomes longer than milliseconds. Some time ago 
Heimburg found that the slow, dominating parts of heat 
capacity and volume-expansion coefficient of DMPC as 
a function of temperature can be superimposed close to 
the melting temperature Tm [H (see also Rcfs. 2 and 3). 
Thus the response functions are related in such a way 
that a single function describes the temperature depen- 
dence of both. 

The fluctuation-dissipation (FD) theorem connects 
(linear) response functions to equilibrium fluctuations. 
The isobaric heat capacity Cp can be calculated from cn- 

] thalpy fluctuations as follows: Cp — {{AH)'^)/{VkBT^), 

< where (. . .) is an average in the constant temperature 
and pressure ensemble and A is deviation from the aver- 
age value. Similarly, volume fluctuations are connected 
to the isothermal volume compressibility by the expres- 
sion kt = {{AV)'^)/{VkBT). If the response functions 
were described by a single parameter , fluctuations are 
also described by a single parameter [l|, 0| and the mi- 
crostates were connected via the relation AH, = j^^^AVi. 
At constant pressure this relation applies if and only if 
AEi = j'"°''AVi (where E is energy), which is the relation 
investigated below. This situation is referred to as a that 
of a single-parameter description . A single-parameter 
description applies to a good approximation for several 



models of van der Waals bonded liquids as well as for 
experimental super-critical argon Q. 

Unfortunately, molecular-dynamics simulations are 
not possible for investigating "single parameter" -ness of 
membranes close to T^, because the relaxation time for 
the collective modes by far exceeds possible simulation 
times. We show below, however, that a single-parameter 
description applies to a good approximation for the slow 
degrees of freedom of the fluid Lq phase, a description 
that applies better upon approaching T„i. At the end of 
this note we briefly discuss how this property may be 
tested in experiments monitoring frequency-dependent 
thermoviscoelastic response functions. 

It is not a priori obvious that a single parameter may 
be sufficient for describing slow thermodynamic fluctua- 
tions of a membrane. For instance, simulations of water 
and methanol showed no "single parameter" -ness. Ap- 
parently, what happens here is that contributions to vol- 
ume and energy fluctuations from hard-core repulsion 
compete with those from directional hydrogen bonds to 
destroy any signiflcant correlation j^. Membranes are 
complex anisotropic systems, and we cannot give any 
obvious reason that volume and energy should correlate 
strongly in their fluctuations. 

A membrane may be pertubated via three thermody- 
namic energy bonds. The change of enthalpy dH can be 
written as a sum of contributions from a thermal energy 
bond, a mechanical volume energy bond, and a mechan- 
ical area energy bond, dH = dE + pdV + HdA, where 
p is pressure, V volume, 11 membrane surface tension, 
and A membrane area. The natural ensemble to con- 
sider is the constant T, p, and 11 ensemble, since the sur- 
rounding water acts as a reservoir. If a single parameter 
controls the microstates, for all states i one would have 
AE, = j'"°^AVi = ^"''^"AAi where the 7's are constants. 
In general, the microstates may of course be controlled 
by several parameters. An interesting question is how 
many parameters are sufficient to describe the membrane 
thermodynamics to a good approximation. This question 
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is addressed below by investigating molecular-dynamics 
simulations of different phospholipid membranes. 

An overview of the simulated systems is found in table 
m The simulated systems include different head groups 
(both charged and zwitterions) and temperatures. All 
simulations was carried out in the constant pressure, tem- 
perature ensemble. The membranes are fully hydrated 
and in the fluid phase. The simulations were per- 
formed using the program NAMD 0] and a modified 
version of CHARMM27 all hydrogen parameter set 
More simulation details are given in Ref. 6. 

The correlations between equilibrium time-averaged 
fluctuations of volume and energy on the nanosecond 
time-scale of a DMPC membrane at 310 K are shown 
on Fig. [TJ If E{t) and V{t) is the energy and volume 
averaged over 1 nanosecond, the figure shows that to a 
good approximation one has 

A;b(I) ~ 7™'Ay(i) (1) 

where 7"°' = aE/(yv is a constant (standard de- 
viation cr) and A is difference from the thermody- 
namical average value. Similar results were found 
for the other membranes studied. Table [J shows 
that volume-energy correlation coefficients {Rev — 

{AEAV)/J{{AEy^) ((AF)2) ) range between 0.81 and 
0.89. 

The correlation depends on the time scales con- 
sidered. This can be investigated by evaluating 
r(t) - CEv{t)/^/ C^{t)Cvv{t) wh ere CAsit) = 
{AA{t)AB{t + t))/^{{AAY){{ABY) is a time corre- 
lation function. T{t) — Q implies that energy at time r 
is uncorrelated with volume at time t -\- t, whereas T{t) 
close to unity implies strong correlation. T{t) is plotted 
in the inset of Fig. [2] At short time (picoseconds) F is 
around 0.5 but approaches unity at t ~ 1 ns. 

The "single-parameter" -ness between volume and en- 
ergy is closely connected to the experimental findings 
of Heimburg since the slow (collective) degrees of 
freedom fluctuating on time scales larger than 1 ns are 
those giving rise to the dramatic changes of the response 
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0.885 
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114 


DMPC 


330 


0.806 


59.0 
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87 


DMPS-Na 


340 


0.835 


45.0 


22 


80 


DMPSH 


340 


0.826 


45.0 


40 


77 



TABLE I: Data from simulations of fully hydrated phos- 
pholipid membranes of l,2-Dimyristoyl-sn-Glycero-3- 
Phosphocholine (DMPC), l,2-Dimyristoyl-sn-Glycero-3- 
Phospho-L-Serine with sodium as counter ion (DMPS-Na) 
and hydrated DMPS (DMPSH). The columns list temper- 
ature, volume-energy correlation coefficient, average lateral 
area per lipid, simulation time in equilibrium (used in the 
data analysis), and total simulation time. 
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FIG. 1: Normalized fluctuations of energy (x) and volume (o) 
for a DMPC membrane at 310 K. Each data point represents 
a 1 ns average. Energy and volume correlate with correlation 
coefficient Rev = 0.89. 

functions approaching T^. Figure [2] shows the time- 
correlation functions of energy, volume, and area of a 
DMPC membrane at 330 K and 310 K. Time constant 
as well as magnitude of the slow fluctuations increase 
when temperature decreases and the phase transition is 
approached. 7''°' in Eq. (P) is 9.3 x 10"^ cm^/J- This 
is of the same order of magnitude as 7""°' — 7.7 x 10^^ 
cm'^/J calculated from the experimental data of Cp{T) 
and K5P'(r) at fl. 

Both volume and energy time-correlation functions 
show a two-step relaxation at 310 K for DMPC (Fig. 
As temperature is lowered towards Tm, the separation is 
expected to become more signiflcant. It makes sense to 
divide the dynamics into two separated processes, a fast 
and a slow collective process. Our simulations suggest 
that the slow degrees of freedom can be described by a 
single parameter, but not the fast degrees of freedom. 

To see the "single-parameter" -ness of the La phase, 
the fast degrees of freedom must be flltered out. Ex- 
periments deal with macroscopic samples where fluctu- 
ations are small and difficult to measure (the relative 
magnitude of fluctuations goes as 1/VN where N is the 



respectively, by definition 
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FIG. 2: Time-correlation functions for a DMPC membrane of 
potential energy Cee (x), total volume Cvv (o), membrane 
area Caa (A), and cross correlation between energy and vol- 
ume Cev {<>)■ Time correlation for membrane areas are scaled 
by a factor 0.2. Panel A shows data at 330 K, panel B at 310 
K. The inset (C) displays r(t) = CEv{t) / y/CEE{t)Cvv{t) 
at 310 K (o) and 330 K (x). F approaches unity at t ~ 1 
ns showing that volume and energy correlate strongly on this 
timescale. 



number of molecules). It is therefore difficult in exper- 
iment to perform the same analysis as we have done 
here; it is easier to measure response functions. Fast 
fluctuations can be filtered out by measuring frequency- 
dependent response functions. The slow collective de- 
grees of freedom give rise to a separate "loss" peaks in 
the imaginary parts. A frequency-dependent Prigogine- 
Defay ratio Axpiuj) was recently suggested as a test quan- 
tity for single-parameter-ness in a paper focusing on the 
properties of glass-forming liquids If Cp{uj), k'^{lo), 
and ap{uj) are the imaginary parts of the frequency- 
dependent isobaric specific heat (per volume), isother- 
mal compressibility, and isobaric expansion coefficient, 



In general Atp(w) > 1, and Atp{lo) = 1 if and only if 
a single parameter describes the fluctuations [4]. The 
quantity 1 / y/T^r^ may be interpreted as a correlations 
coefficient. 

In summary, we found strong volume-energy corre- 
lations of the slow degrees of freedom in molecular- 
dynamics simulations of different phospholipid mem- 
branes in the phase. An experimental test was sug- 
gested. 



ACKNOWLEDGMENTS 

The authors would like to thank Nick P. Bailey, 
Thomas Heimburg, and S0ren Toxvaerd for fruitful dis- 
cussions and useful comments. This work was supported 
by the Danish National Research Foundation Centre for 
Viscous Liquid Dynamics "Glass and Time." GHP ac- 
knowledges financial support from the Danish National 
Research Foundation via a grant to the MEMPHYS- 
Center for Biomembrane Physics. Simulations were per- 
formed at the Danish Center for Scientific Computing at 
the University of Southern Denmark. 



[1] 
[2] 



[3] 
[4] 

[5] 
[6] 

[7] 



T. Heimburg, Biochim. Biophys. Acta 1415 (1998) 147. 
T. Heimburg and A. D. Jackson, Proc. Natl. Acad. Sci. 
102 (2005) 9790; T. Heimburg, Biochim. Biophys. Acta 
1415 (1998) 147. 

H. Ebel, P. Grabitz, and T. Heimburg, J. Phys. Chem. B 
105 (2001) 7353. 

N. L. EUegaard, T. Christensen, P. V. Christiansen, N. B. 
Olsen, U. R. Pedersen, T. B. Schr0der, and J. C. Dyre, J. 
Chem. Phys. 126 (2007) 074502. 

U. R. Pedersen, N. P. Bailey, T. B. Schr0der, and J. C. 
Dyre, ArXiv: cond-mat/0702146 (2007). 
U. R. Pedersen, C. Leidy, P. Westh, and G. H. Peters, 
Biochim. Biophys. Acta 1758 (2006) 573; U. R. Pedersen, 
G. H. Peters, and P. Westh, Biophysical Chemistry 125 
(2007) 104. 

J. C. Phillips, R. Braun, W. Wang, J. Gumbart, E. 
Tajkhorshid, E. Villa, C. Chipot, R. D. Skeel, L. Kale, 
and K. Schulten, J. Comput. Chem. 26 (2005) 1781. 
N. Foloppe and A. D. MacKereU, J. Comput. Chem. 21 
(2000) 86. 



